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Abstract. We investigate the behavior of the thermoelectric power 5* in disordered systems 
close to the Anderson-type metal-insulator transition (MIT) at low temperatures. In the 
literature, we find contradictory results for S. It is either argued to diverge or to remain 
a constant as the MIT is approached. To resolve this dilemma, we calculate the number 
density of electrons at the MIT in disordered systems using an averaged density of states 
obtained by diagonalizing the three-dimensional Anderson model of localization. From the 
number density we obtain the temperature dependence of the chemical potential necessary 
to solve for S. Without any additional approximation, we use the Chester-Thellung-Kubo- 
Greenwood formulation and numerically obtain the behavior of S at low T as the Anderson 
transition is approached from the metallic side. We show that indeed S does not diverge. 
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1 Introduction 

In this paper, we study the low temperature behavior of the thermoelectric power S in 
disordered systems near the Anderson-type metal-insulator transition (MIT). In the 
framework of linear response theory, S, commonly abbreviated as the thermopower, 
is the coefficient that relates the temperature gradient in an open circuit with the 
induced electric field. In the metallic regime, the Sommerfeld theory states that S 
is directly proportional to the negative temperature —T But at a disordered- 
induced MIT, such as the Anderson transition in three dimensions (3D) [||, it is still 
not a settled issue how S behaves at low T. Theoretical studies have either claimed 
that it diverges or that it remains a constant as the MIT is approached from 
the metallic side at low T. Moreover, comparing the results of the latter theory 
with that of experiments conducted on doped semiconductors |^ and on amorphous 
alloys 1^ shows that measurements of S are two orders of magnitude higher than 
those predicted in theory. Thus, it is of great interest to investigate the behavior of 
S at low T near the Anderson-type MIT. Here, for simplicity, we consider only the 
diffusion part of S, that is, we consider only the electronic contribution and neglect 
any possible electron-phonon interactions. In addition to S, we shall also compute 
thermal transport properties such as the thermal conductivity K and the Lorenz 
number Lq. 
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2 Theoretical background 

The derivation of the thermopower is based on the kinetic coefhcients of the Chester- 
Thellung-Kubo-Greenwood formulation of linear response 0, 

= {-ir^ p^ A(E)[E ^ ^,{T)X+^-^ dE ^,J = 1,2 (1) 

where E is the energy, A{E) contains all the system-dependent features, is the 

chemical potential, f{E,ii,T) = 1/[1 + exp{[E — fi{T)]/kBT)] is the Fermi function, 
and ks is the Boltzmann constant. These coefficients relate the electric field e, con- 
centration gradient V/z and temperature gradient VT to the expectation values of 
the induced electric (ji) and heat (j2) current densities 

(j^) = [~La{y^i ~ W) ~ LaT'^VT\ . (2) 

where e is the electron charge. Measured under the assumption that there is no 
electric current and concentration gradient, the thermopower is thus given as 

\e\TLii 

The Anderson transition is then incorporated into the measurement of S by setting 
the function A{E) in the coefficient Lij as proportional to the critical behavior of the 
d.c. conductivity a at the MIT, that is, 

r a\E-Ecr E>Ec 
"^^^^-I E<Ec 

where a is a constant, v is the conductivity index and Ec is the mobility edge. 
With this assignment the coefficient Ln is simply a. Furthermore, since K is the 
coefficient that relates the temperature gradient to the induced heat current, it's 
low temperature behavior at the MIT can be determined in a similar manner from 
with the assumption that there are no particle currents, and using the Anderson 
transition form of A(E) as given above. Then the Lorenz number Lg = {e/kB)'^<j/KT 
quickly follows. Thus, the low T behavior of 5", K and Lq at the Anderson transition 
follows easily after obtaining the kinetic coefficients, Eq. (|l|). 

2.1 Divergent thermopower 

A divergent S at the Anderson transition E = Ec is obtained if one uses the Sommer- 
feld expansion to get the low-T leading contribution to Lij Q . This method assumes 
that the chemical potential fi is equal to the Fermi energy Ep even for finite T. How- 
ever, /i = Ep only at T = |jl|]. A more serious approximation of the Sommerfeld 
expansion is the assumption that A{E) is a smoothly varying function ai E = Ec- 
This is not the case at the Anderson transition, as can be readily seen in Eq. (|^). 

2.2 Fixed-point thermopower 

The approach proposed by Enderby and Barnes evaluates the kinetic coefficients 
at /J, = Ec for finite T, and afterwards the limit T ^ is taken. They find that the 
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thermopower is a constant at the mobility edge for T — > 0, and is given by 

^^_k_s{j^I^ (5) 

|e| V ly 

where h = In 2, = (1 - 2(i-'^))r(i/)C(i^) for ReH > 0, / 1, with T{v) and C{v) 
the usual gamma and Riemman zeta functions. Hence, S solely depends on v. 



3 Calculation of the temperature dependent thermopower 

One can determine the temperature dependence of the thermopower if one knows 
how /i varies with T . This information can be obtained from the number density n of 
electrons at the MIT. In general, for any set of noninteracting electrons, the number 
density is defined as 

/oo 
dEg{E)f{E,fi,T) , (6) 
-C30 

where g{E) is the density of energy levels per unit volume. Using the above equation, 
we numerically calculate n using an averaged density of states g{E) obtained by 
diagonalizing the Anderson model of localization. Earlier, we determined the averaged 
density of states for a 3D isotropic Anderson model with disorder = 12 Note 
that since our objective is to compare our theoretical results for S with experimental 
measurements, such as those from amorphous alloys, the hopping parameter t is of 
the order of 1 eV. Hence, we have expressed all energy units in terms of t unless 
otherwise specified. We have selected the value of W to be strong enough, such that 
we do not have singularities in the density of states. Yet, it should not be too strong, 
i.e. too close to the critical disorder. For this particular value of W, the value of Ec is 
approximately —7.5, according to the mobility edge trajectory Ec{W) calculated in 
Ref. The conductivity index j/ is « 1.3, according to a current numerical estimate 
[|o| . Then we integrate the density of states for E < Ep to obtain the corresponding 
value of n for a given value oi Ep s^t T — 0. Keeping n fixed at this value, we vary 
T in Eq. (|^) and numerically determine the variation of /i. Using this information in 
Eq. (|l|) , we solve for Lij . It is then straightforward to determine S for a particular 
value of Ec from Eq. (H) . 



4 Results and discussion 

In Fig.|l[ the temperature dependendence of the chemical potential is shown together 
with the averaged density of states from which it was measured. Note that from this 
smooth density of states, we obtain a T dependence of fi which barely changes when 
one selects Ep in the metallic or the localized region. However, its slope changes 
much faster as compared to the chemical potential from a free electron gas as shown 
in Fig. 1. Note that this free electron result was also similarly obtained from the same 
expression for n given in Eq. (^, but using the Sommerfeld expansion in order to 
obtain /x. 

Next, Fig.|^ shows our thermopower measurements. The curves at the top of Fig.|^ 
clearly show the MIT, the dividing line between the metallic {Ep > Ec) and localized 
{Ep < Ec) regions. As T ^ 0, S' gets more negative in the localized region, while 
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Fig. 1: Top: The low T behavior of ^. Near the MIT, /i(T) is similar in both the localized 
and the metallic regions. Bottom: The averaged density of states of a 3D isotropic Anderson 
model with W — 12. For clarity only every 10th data point is marked by a symbol (o). 



5 — > in the metallic region. As we move further away from the MIT towards the 
metallic region at low T, S behaves as expected from the Sommcrfeld theory, that is, 
linearly proportional to —T. This indicate nonzero values of a{Ep) confirming the 
metallic nature in this energy region. More importantly, we see that S* is a constant 
at the MIT, Ep = Eq. As T ^ 0, it approaches the value -228.4 ^iV/K. This value 
agrees with the T-independent value for v — 1.3 as predicted by Eq. (||). At the MIT, 
a negative S value of the order of hundreds of fiV/K has never been experimentally 
observed to the best of our knowledge. To see the T-independence of S at the MIT, 
we refer to the bottom of Fig. ||. Here we show the behavior of S at different Fermi 
energies for different temperatures. It is clearly demonstrated in the inset that for 
different values of T, S is a. fixed point at the MIT (/i = Ec) verifying what Enderby 
and Barnes had previously concluded [Q. 

Similarly, we have studied the other thermal transport properties K and a. Our 
preliminary investigation shows that K —> as T —> at any energy region. Fur- 
thermore, in the metallic phase, Lq approaches the value 7r^/3 which according to 
the law of Wiedemann and Franz is a universal value for all metals (see for example 
Refs. [0, 0]). At the MIT, however, Lq has a value dependent only on the conductivity 
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Fig. 2: Top; The low T behavior of S. Note that 5" docs not diverge at the MIT as T 0. 
Bottom: Same data as in the top panel, plotted with respect to /i for different Fermi energies. 
The lines connect isotherms of S. As shown in the inset, S is a fixed point at the MIT. 



index v. Detailed results of these transport properties will be discussed elsewhere. 



5 Conclusions 

In this work we have studied the low temperature behavior of the thermoelectric 
power for the 3D isotropic Anderson model close to the MIT. We have numerically 
obtained the temperature dependence of the chemical potential necessary to solve for 
5* from the general expression of the number density for any set of nonintcracting 
electrons. We have shown that is very similar regardless which energy region 

close to the MIT one considers. Using this result and the Chester-Thellung-Kubo- 
Greenwood formulation, our calculations yield a sharp contrast of the S behavior 
between metallic and localized regions clearly outlining the MIT. Finally, as the MIT 
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is approached from the metahic side 5 is a fixed point. As T ^ at the MIT, S 
approaches the fixed-point value predicted by Enderby and Barnes which for u = 1.3 
is S = —228.4 fiV/K. Therefore, we have estabhshcd that as the MIT is approached 
at low T the thermopower does not diverge but remains a constant. Its fixed-point 
value depends only on the critical behavior of a. How S behaves for varying degrees 
of disorder is a subject of further investigation. 

We thank T. Vojta for helpful (iiseiissions. We also gratefully acknowledge financial support 
by the DFG through Sonderforschungsbereich 393. 
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